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We investigate the interaction between a single mode light field and an elongated cigar shaped 
Bose-Einstein condensate (BEC), subject to a temporal modulation of the trap frequency in the tight 
confinement direction. Under appropriate conditions, the longitudinal sound like waves (Faraday 
waves) in the direction of weak confinement acts as a dynamic diffraction grating for the incident 
light field analogous to the acousto-optic effect in classical optics. The change in the refractive index 
due to the periodic modulation of the BEC density is responsible for the acousto-optic effect. The 
dynamics is characterised by Bragg scattering of light fom the matter wave Faraday grating and 
simultaneous Bragg scattering of the condensate atoms from the optical grating formed due to the 
interference between the incident light and the diffracted light fields. Varying the intensity of the 
incident laser beam we observe the transition from the acousto-optic effect regime to the atomic 
Bragg scattering regime, where Rabi oscillations between two momentum levels of the atoms are 
observed. We show that the acousto-optic effect is reduced as the atomic interaction is increased. 

PACS numbers: 03.75.Kk, 05.45.-a, 32.80.-t, 42.50.-p, 47.54.-r 

I. INTRODUCTION 

When a liquid traversed by compression waves of shorter wavelength is irradiated by visible light, diffraction 
phenomenon similar to that due to a grating is produced. This phenomenon waspredicted by Brillouin [l[ in 1921 
and subsequently verified experimentally by Debye et al. 0, and Lucas et al. This effect is considered to be 
an important discovery since this led to the interesting application of acoustic modulation of optical radiation by 
naturally birefringent crystal. The linear acousto-optic (also known as photo-elastic) effect involves the first order 
changes in the optical properties of insulators due to acoustic strain. The acousto-optic effect, however has not yet 
been suggested for the most recently created state of matter, the Bose-Einstein condensate (BEC). 

An important question in this context is the generation of compression waves in a BEC which can act as a dynamic 
diffraction grating for the incident light held. The generation of such a kind of longitudinal sound like waves was first 
predicted by Staliunas et al. 0, Q and experimentally achieved recently in an elongated BEC in the direction of weak 
confinement by periodically modulating the radial trap frequency Q . The waves thus generated are known as Faraday 
waves. The modulation of the radial trap frequency leads to a periodic change of the density of the cloud in time, 
which is equivalent to a change of the nonlinear two-body interaction. This in turn, lead to the parametric excitation 
of longitudinal soundlike waves in the direction of weak confinement. In this work, we report on the diffraction of 
light waves from a BEC subjected to a temporal modulation of the trap frequency in the tight confinement direction. 
In contrast to the classical acousto-optic effect in liquids and crystals, the two-body atom-atom interactions and the 
back action of the light fields on the atomic medium in the present system generates new and interesting physics. 
The light fields create an effective optical grating inside the BEC which Bragg scatters the very atoms which are 
generating the matter-wave grating for the light field. This complex interplay between the atoms and the photons 
gives the acousto-optic effect a quantum character unseen in the classical case. The atomic interactions are tunable 
by Feshbach resonances and this allows us also to investigate the effects of the mean-field atomic interaction on the 
diffraction process. 

In the context of interaction of light fields with BEC, the long coherence time of the BEC offers the possibility 
to study collective optical effects of atoms in atom optics such as Bragg scattering of atoms from standing wave 
laser field and superradiant Rayleigh scattering. The photon exchange between different ultracold atoms result in a 
long-range inter-atom interaction which under appropriate conditions leads to an effective "Kerr- type" nonlinearity 
of atomic waves 0. Analogous to nonlinear quantum optics, multiwave mixing has been predicted Q and observed 
[9l. fTod in BEC. A nonlinear quantum optics approach to the theory of optically driven BEC's was given by [ll|, where 
it was shown that quantum noise can generate entangled atom-photon pair which exhibit nonclassical correlations 
similar to those seen betwe en p hotons in optical parametric amplifier. The recent observation of superradiant Rayleigh 
scattering from BEC's 12, 1 1 31 ] has significantly extended our knowledge about collective emission processes. The basic 
theory of superradiant Rayleigh scattering has been discussed by a number of authors p^-[l7|. 
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Figure 1: Schematic diagram showing the Bragg scattering of light incident at an angle 9i with wavevector ki by Faraday waves 
of wavevector k s in an elongated cigar shaped BEC. The incident light is diffracted at an angle 8d and a wavevector kd- The 
polarization of the incident light is chosen to be in the plane defined by the z axis and the direction of the wave-vector of the 
incident light so that the superradiant Rayleigh scattering is suppressed. 



II. FARADAY PATTERN AS A DYNAMIC DIFFRACTION GRATING 



In our theoretical treatment, we consider an elongated cigar shaped condensate consisting of N two-level atoms, 
oriented along the z axis. The trapping potential of the ground state is periodically modulated and is written as 

I [J r 1 _ 

Vi(r, t) = — (*) (x 2 + y 2 ) + Q 2 z 2 , with fij_ = £l± (1 + acosQt). Here f2j_ and Q|| are the trap frequencies along 

the radial and the longitudinal (z) direction respectively. Also, Ct± is the amplitude of the radial trap frequency and 
Q is the frequency of the trap modulation. The condition f2j_ » f2|| is assumed. Due to the periodic modulation of 
the trap, the ground state wavefunction is also modulated as: 



ip 1 {r,t) = ip (r,t){l + w(t) cos k s .fj. (1) 

Here, %po(r,t) is the wavefunction of the homogeneous BEC in the absence of any modulations and w(t) is the 
complex valued amplitude of the perturbation. The periodically modulated BEC is exposed to a far off-resonant laser 
field whose polarization is parallel to the longitudinal axis (z). Under proper conditions, part of the input optical 
beam is diffracted into a new direction while simultaneously being shifted in frequency and wavenumber by an amount 
equal to the frequency of the Faraday wave (lo s ) and wavenumber (k s ) respectively. Unlike the classical acousto-optic 
effect, here the choice of the polarization of the incident light is extremely important since the dominant competing 
process superradiant Rayleigh scattering can completely mask the acouto-optic effect in a BEC. According to ref. [l2l ]. 

dN- sin 2 

the growth rate of Nj recoiling atoms is — — ^ = RNq —^-H,j(Nj + 1), where R is the rate for single-atom Rayleigh 

at 87T/3 

scattering, which is proportional to the laser intensity, iVo is the initial number of atoms, 9j is the angle between the 
polarization of the incident light and the direction of emission (z axis) and £lj is the solid angle. For polarization 
parallel to the z axis (9j = 0), light emission into the end- fire mode is suppressed and the higher order momentum 
side modes of the condensate are not generated and the distribution of the atoms follow the dipolar pattern of normal 
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Rayleigh scattering. For perpendicular polarization (6j = 7r/2), photons are emitted along the z axis and the recoiling 
atoms appear as highly directional beams at an angle 45° with respect to the z axiz. This is the superradiant Rayleigh 
scattering which can destroy the acousto-optic effect which we wish to study in this work. Due to the highly directional 
nature of the momentum transfer in the superradiant scattering event, the kinetic energy gained by the atoms due to 
the recoil , far exceeds the kinetic energy gained due to the periodic modulation of the trap. In the case of normal 
Rayleigh scattering, due to the random nature of the spontaneous emission, the direction of the recoil is random, 
leading to momentum diffusion. The amount of atoms lost from the original condensate in the superradiant Rayleigh 
scattering is extremely high compared to that from normal Rayleigh scattering fl2j. In order to suppress super-radiant 
Rayleigh scattering, the incident light beam is taken to be linearly polarized in the plane defined by the condensate 
axis and the wave vector of the light [HI, H3] • Note that experimentally for parallel polarization and small detuning, 
yet another new effect is produced namely the superradiant Raman scattering [Til . In the superradiant Raman 
scattering, only two modes of the recoiling condensate atoms are populated and the decoherence is also much faster 
compared to superradiant Rayleigh scattering due to inhomogeniety in the magnetic field. This fast decoherence rate 
can be traced to a new type of coherence grating that is formed and this grating is highly sensitive to magnetic fields. 
A slightly inhomogeneous magnetic field could be used to suppress the coherence grating which could also interfer 
with the quantum acousto-optic effect. There is yet another effect which is in phase with the acousto-optic effect, 
the Bragg scattering of the condensate atoms from the o ptic al lattice created due to the interference between the z 
component of the incident and the diffracted light waves [21(. In section IIB, we treat the light field classically with 
fixed intensity and treat the atoms in the mean field regime. This treatment is useful to gain insight into the dynamics 
of the BEC density and in particular the time evolution of the shape of the cloud. In section III, we take the light as 
quantum and this treatment helps us to study the time evolution of the light fields. We will show in section III that 
the Bragg scattering of the atoms can be minimized if the decay rate of the light fields are greater than the atom-field 
coupling strength (this can be controlled by the intensity of the light field). 

A. Bragg scattering condition 

Figure 1, shows the schematic setup with which one can observe the acousto-optic effect in an elongated cigar shaped 
BEC. The Faraday waves act as dynamic diffraction grating and Bragg scatter the incident light. Let the incident 
and the diffracted optical field consist of two plane waves with frequencies u>i and u>d respectively. The two waves 
propagate along ki and k c i.e, respectively, where ki in general is not parallel to kd,i- There are two main diffraction 
regimes, the Raman-Nath regime and the Bragg regime (22j . When the parameter x — ekL/k s sin9i > 1, we are in 
the Raman-Nath regime which is characterized by the presence of many diffraction orders and when \ < 1, we are 
in the Bragg regime which is characterized by the presence of only the first order diffraction. Here e and fcj, are the 
dielectric constant of the BEC and the wavenumber of the laser beam. Energy and momentum conservation leads to 
the following relations: 

Ud,e —Ui = ho a (2) 

kd,e — ki = £k s (3) 

Here I = 0, ±1, ±2... is the order of diffraction and kd.t is the diffracted wavevector of the I th order. The interaction 
can be viewed in the following terms. A photon with energy huii and momentum hki is incident on a Faraday wave 
of frequency uj s and wave momentum hk s . The incident photon and a phonon of the Faraday wave are annihilated, 
giving rise instead to a new photon at LOdi , kd t e- 

The Eqn. [3] yields the usual familiar diffraction condition for fcj = kd,e = k^: 

k L (sin ~ sin &i ) = ik s (4) 

The periodic nature of the Faraday waves play a role similar to that of a regular arrangement of atomic planes. In 
Eq.(4), we have assumed fc; = kd = ki,. Also fcj, 10 6 m -1 and k s « 10 5 m _1 . 

B. BEC dynamics in the presence of classical light field and Faraday instability 

The system consisting of the periodically modulated condensate with N two-level atoms and the light fields is 
described by the Hamiltonian 
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h 2 V 2 

2m 



+ V t (r, t) ) i)j (f, t) + I d 3 ri>l (r, t)hu a ^ 2 (r, t) + H L 



(5) 



Here, Hl is the Hamiltonian for the lasers, w a is the atomic transition frequency and /z is the chemical potential. 
f2j and fi^ are the Rabi frequencies of the incident and the £th order diffracted laser beams. The £th order Rabi 

frequency is related to the incident Rabi frequency as fld,i=^i exp{i£ir /2) Je(ekLL/2\/ n 2 — sin 2 Oi) [2^ |. Je is the £th 

order Bessel function , n is the refractive index and L = 10/im is the thickness of the condensate. U = is 

m 

the two body interaction parameter, where a s is the interatomic s-wave scattering length (taken to be positive) and 
m is the mass of the atom. 

We assume that the population in the excited state is extremely small so that the atom-atom interaction in the 
excited state is neglected. Here, ipi and ip2 are the ground and excited state field operators respectively. ki(u>i) and 
kd(oJd)is the wavenumber (frequency) of the incident and diffracted light respectively. V$ (i = 1, 2) is the trap potential 
of the ground state (i — 1) and the excited state (i = 2). 

We can now write down the Heisenberg equation of motion for atomic field operators, and ip2 and then eliminate 
the excited state under the condition that the detuning 6 = LUi.d — u a is much larger than the natural line width of 
the atomic transition. The kinetic energy and the trap potential of the excited state is dropped under the assumption 
that the lifetime of the excited atom, which is of the order 1/8, is so small that atomic center of mass motion may be 
safely neglected during this period. This yields the following equation of motion (the Gross-Pitaevskii equation) for 
the ground state in the mean field regime where we replace the atomic field operators by the corresponding c numbers: 



ih 



dt 



2m 



+ Vx{r,t) + V ov {r,t) -fi 



(6) 




E 



dj\ 



(7) 



Where, Aw{ = ujd,e — = £uj s and Akf — kd.i — ki = £k s are given by the Bragg diffraction condition [21 Here 
V\ — m/2(fl 2 L (t)(x 2 + y 2 ) + fl^z 2 ). The first two terms in Eqn.([7]) is simply the static optical potential while the 
next two terms are the optical potential formed due to the interference between the z component of the incident and 
diffracted light. The self consistent optical potential thus formed has the same spatiotemporal variation as the Faraday 
wave. Actually, if the decay rate of the light field is much more than the parameter QiQ,d,e/$ then the optical potential 
formed will have a negligible influence on the acousto-optic effect. The chemical potential /i in the Thomas- Fermi 
approximation is determined by the number of particles, the scattering length a s and the trap parameters as [24| 
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hio. 



(8) 



Where uj = (ujxLOyLOz) 1 ^ 3 , a mean = y/h/mu) , u> x ,u> y and lu z are the trap frequencies along the x,y,z directions 
respectively and a s is the scattering length for Rubidium atoms. All lengths are scaled with respect to a± = y/h/mfl± 
and all frequencies are scaled with respect to f2j_, where Oj_ is the amplitude of the trap frequency along the radial 
direction and m is the mass of the atom. The spatiotemporal oscillations of the BEC density along the weakly 
confined space occur at half the trap modulation frequency (oj s = 0/2) . Here we are neglecting the damping of the 
Faraday waves. This can be approximately true if we have a strong and constant driving so that a steady Faraday 
pattern is maintained during the process of diffraction. Inclusion of damping only changes the threshold for pattern 
formation. We take the experimental conditions described in Ref. [6(. Specifically, we take the number of condensed 
87 Rb atoms as N = 5 x 10 5 contained in a magnetic trap with frequencies {f2j_/(27r), f2ii/(27r)} = {160.5, 7} Hz where 



5 



o 



> 

T3 



05 

s 
Q 



0.12 
0.10 
0.08 
0.06 
0.04 
0.02 
0.00 




10 20 30 40 50 60 
r (dimensionless time) 



Figure 2: Plot of deviation from initial density for £ — 0, ±1, p, = 2, a; = 2, a = 0.2 and for two values of the strength of the 
optical potential f2 | 2 /2<5f2x = 0(thick line) and |f2°| 2 /2<5f2x = 0.8(thin line). The solutions are found to be relatively more 
unstable in the presence of the light field. 

the radial trap frequency = f2j_ (1 + acosfii) has a modulation of 20 percent (a = 0.2) and a modulation 

frequency = uj = 2. In terms of the scaled variables, Eqn.(6) is reduced to a dimensionless one dimensional 

equation (this can be done by separating the wavefunction into a radial and a longitudinal part Q). The radial part 
is approximately taken to be Gaussian and the radial part is integrated out and we are left with an effectively one 
dimensional equation: 



1 dr 



dZ 2 



Z 2 + V op (Z, T )-fi 



<j) (Z) + ^\cb(Z)\ 2 ( t>(Z) 



(9) 



V op (Z,r) 



I a- 



E 



a, 



-ik d Z 



-iu) s r AkdZ 



(10) 



Here jl = 2(/x— a>j_) and the dimensionless parameters arc defined as : \x — /i/Q,±, lu±_ = f2_i_/fi_i_, = w s /r2j_, = 
f2||/f2j_, r = Sl^i, Z = z/a± and kd/cL±- The above equations are solved numerically and the resulting normalized 
deviation from initial density (scaled central density) is plotted in Fig. 2 for £ — 0, ±1, jl = 2, u> = 2, a = 0.2 and for 
two values of the strength of the optical potential |f2 | 2 /25f2j_ = 0(thick line) and |r2°| 2 /2<5f2^ = 0.8(thin line). Here 
we have taken flj and 0^ j as real and Hi = \/3&ld,e = The solutions are found to be relatively more unstable in 
the presence of the light field since in the presence of the light field, the amplitude of the perturbation grows much 
faster. The faster growth of perturbations in the presence of the light field is probably due to excitation of the several 
other modes which makes the system unstable. Actually for a constant driving the condensate behaves like an impact 
oscillator which is accompanied by increasing transverse cloud radius with time [||. 

III. TWO-LEVEL MODEL OF THE QUANTUM ACOUSTO-OPTIC EFFECT 



In this section, we adopt a semi-classical approach to investigate the problem of Bragg scattering of light from 
Faraday waves travelling along the long axis of an elongated cigar shaped BEC. The Faraday pattern consists of 
Faraday waves travelling in opposite direction with wavevectors k s ,—k s . We consider the situation as shown in Fig. 
1 where an incident photon of frequency u>i and wavevector ki is scattered by a Faraday wave phonon of frequency uj s 
and wavevector k s travelling in the positive z direction. The energy and momentum conservation condition requires 
that the diffracted photon have a frequency u>i +tu s and wavenumber ki + k s . The interaction in the case of scattering 
of incident photons by Faraday wave travelling in the negative z direction can be viewed in the following terms. A 
photon with energy fkOi and momentum hki is incident on a Faraday wave of frequency oj s and momentum —Hk s . 
The incident photon is annihilated, giving rise instead to a new photon at u>d, fc<2 and a phonon u> s , k s travelling in 
the positive z direction. The above two interactions effectively involves only three momentum levels i.e. level with 



6 



zero momentum and that with ±hk s . The loss due to decoherence and decay of the light field is taken into account. 
The above process is described by the quantized Hamiltonian: 

2 f ( h 2 V 2 \ f 

H = / d 3 rtp](f, t) ( — — + Vi(r, t) J ^(f, t) + d 3 ft/4(r, t)hoj a ip 2 (r, t) + hu)i&\ai + frw 2 a\a 2 



-iftjsi&i y d 3 rMr,t)e^l(f,t) +92 a 2 J dV^(r, t)e Jfcd - r Vi(r, t)| + H.c 
+U { d 3 f^\{r,t)^l{r,t)Mr,t)^x{r,t) (11) 



Here d\ and 02 are the destruction operators for the incident and the diffracted light respectively and g\ and 172 are 
the atom field electric-dipole coupling constant for the incident and the diffracted light fields respectively. The ground 
and excited state atomic field operators are tpi and 1P2 respectively. We can write down the equation of motion for 
the operators 4>2, 0,1 and d 2 from the Hamiltonian 1111 

ih ^T = (-^T- + V i) & + ^9id\e-^^ 2 + ihg 2 d\e-^-^ 2 + U\^\ 2 ^. (12) 
ot \ 2m J 

ih ^T = (-^T- + v *) & - - ihg 2 a 2 e i% "-% + hw a ^ 2 . (13) 

at \ 2m J 

ih^p- = huji ai + ihgi f d 3 r^Je~**** r ^2 — (14) 
at j 

= huj 2 d 2 + ihg 2 / d 3 r^|e _4fcci - , 'V2 - nd 2 . (15) 

Here, k ~ c/2L is the phenomological radiation loss term taken to be same for both the incident and the diffracted 
beam and L is the length of the condensate [2(|. Introducing the operators ip'^ = ^ 2 e~' l< - UJi+u ' d ^ 2t and bj = cije~ tuJit , 
j = 1,2, we adiabatically eliminate the excited state. We now replace the bosonic operators in the above equations 
by the coherent wavefunction of the condensate < ipj >= V'j (mean-field approximation) and the classical light fields 
amplitudes < bj >= bj (in the case of large intensity of the light fields). Integrating out the radial dimension, we get 
the following equations of motion for the new variables in one dimension: 

iH ^t = " Ml + ~5 + 92b * 2b2 + ^i^-^V^ + gig^b^e-^] I ^ 1 + Um 1 \ 2 ^ 1 , (16) 

= -Mh - iG'he" 1 ^ [ dze ih ' z \i\) X ? - Kb 1: (17) 
dt o J 

^ = -i4h - iG'he™' 1 [ dze-^ z \^\ 2 - nb 2 , (18) 
at o J 

Where, G = ^j^-ir I Q± ] , m = p, - m± and U{t) = U{1 + acosfit). The Eqns. US Q7| and Ql are the 
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coupled matter-wave and the light fields which determine the complete dynamics of the system. The third term on 
the right-hand side of Eqn. [T|5] represents the optical potential (with the same periodicity as that of the Faraday 
wave) created by the interference between the z component of the incident and the diffracted light fields and whose 
amplitude depends on time according to Eqns. rT71and[T51 This optical grating Bragg scatters the condensate atoms. 
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Figure 3: Time evolution of the normalized population difference(W(t), thin curve) between the two momentum modes and 
the intensity difference (N p (t), thick curve) between the diffracted light beam and the incident beam, (a) The parameters 

IO°|2 

chosen are: g = 10, g' = .1, a = 0.1, k = 0.1 = 0.1, 7 = 0.1, u = 2, N = 10 5 , a s /a mean = 0.00016, -UJ- = 0.1, 

nil 7 

■=- li - = . The population in the two momentum modes undergo a damped Rabi oscillations due to Bragg scattering of the 

Q± 160 

atoms from the optical grating when the light-atom coupling is greater than the loss rate of the light field. In this regime, the 
Rabi oscillations has a small contribution from the quantum acousto-optic effect and this is visible as a small beat frequency in 
the oscillations. The intensity difference, N p (t) also undergoes oscillations and approaches a steady value after some time, (b) 
In the opposite limit when the decay rate of the light field (k = 1.1) is larger compared to the light-atom coupling (g — 1.0), 
the Rabi oscillations are highly damped and the initial population difference decays to a stationary value, (c) On increasing 
the particle number N = 5 x 10 5 , the amplitude of the Rabi oscillations decreases, (d) When the amplitude of the perturbation 
is increased, a = 0.5, the beat frequency in the Rabi oscillations becomes strong due to an increased contribution from the 
quantum acousto-optic effect. 



The second terms on the right-hand side of Eqns. [T7] and [TS] are the matter-wave Faraday grating. An important 
point to note is that in deriving the EqnsITTJl [T71 and [TSl we have assumed the Bragg condition 2J If the condensate 
(~ 125/ito) is much larger than the Faraday wavelength(10/xm) then periodic boundary condition can be assumed on 
k s z and the wavefunction can be written as a Fourier series [261 ] 

Mk s z) = c n {t)u n {k s z)e- m ^ t . (19) 

n 

Where u n {k s z) = —j=exp(ink s z) are the momentum eigenfunctions with eigenvalues^ = n(hk s ). If we work close 

to the Bragg regime and the condensate is homogeneous along the long axis (this is true for a weak confinment along 
the long axis), we can assume that the only two momentum levels involved in the process are n and n + 1 (the two 
level approximation is true only if one diffracted beam is considered) . We get the following equations for the variables 
W = |c„| 2 - |c„+i| 2 , S = c„< +1 , A = ihb*, N p = b*b 2 - b$h 
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dS_ 

dr 



^ = -2g(A*S + AS*) (20) 



-i(A n -U'(T)W)S + gAW--fS (21) 



dA 

— = g'SN p -kA + in' A (22) 
dN v g' dW 

-W = (23) 

In the above equations, all the physical quantities have been made dimensionless according to the scaling given in 
the previous section. Here, g = 2 , g' — G'/fij_,A„ = (S n — 5 n+ i)/£l±, S n = — — nu; s , U'(t) = UM/MlxV 

g 2 - g 2 

(V is the volume of the condensate), k' = — r= and n — > k/Q±. We have introduced by hand 7 (dimensionless 

w.r.t n±) as the decay rate of the atomic coherence between the two motional states n and n + 1. The decoherence 
is induced by spontaneous emission, phase diffusion, Doppler and inhomogeneous broadening [26|. The dynamics 
of the system shows a phase matching between two kinds of phenomena, the acousto-optic effect and the Bragg 
scattering of atoms from the optical potential formed due to the interference of the incident and the diffracted light. 
The usual acousto-optic effect is dominant when either the atomic decoherence rate 7 or the radiation loss term k 
are greater than g and g' and the Bragg scattering of the atoms is dominant mechanism in the opposite limit. In the 
absence of interactions (U = 0), A„ = is the atomic Bragg scattering condition, arising from momentum and energy 
conservation. The two motional states assumed to be involved in the acousto-optic effect are the n = (stationary 
condensate state without any Faraday wave) and the m = n + 1 = 1 (state with one unit of Hk s ). 

The above equations are solved numerically in FigJ3£i for conditions g = 10, g' = .1, a = 0.1, k = 0.1 ,k' = 0.1, 

|fi°l 2 fin 7 

7 = 0.1, lu = 2. N = 10 5 , a s /a mea n = 0.00016, ^-r^- = 0.5, = . The initial conditions used are W(0) = 

ou± f2j_ 160 

0,iVp(0) = — 1,A(0) — 0, 5(0) = 0.1. We found that the conclusions that follow does not depend on the initial 
conditions, however a small but finite value of S(0) is required to observe the dynamics (this implies that an initial 
coherence between the momentum modes is needed for the dynamics to start ). The diffracted light field grows in 
time and reaches a steady value and follows W(t). As the diffracted field intensity grows, the population of the lower 
momentum state |c„| 2 also grows since the condensate in the higher momentum state relaxes to the lower level after 
transferring momentum to the light field. After some time due to the Bragg scattering, the atoms perform a weakly 
damped Rabi oscillation caused by Bragg transition between the two momentum states p — and p = Hk s . This 
process continues till a steady state is reached by both the atomic population difference (W(t)) and the intensity 
difference N p (t), due to a finite decay of the atomic coherence and the light field. Note that till time t = 60ms, only 
about less than half the energy is transferred to the diffracted beam due to the Bragg scattering of atoms. The rate 
at which atoms are transferred between the two momentum modes due to the Bragg scattering is different from that 
due to the acousto-optic effect and this results in the presence of beats in the Rabi oscillations in FigJSk.. 

The Rabi oscillations are suppressed if the decoherence factor 7 or the light field decay rate k is greater than g. 
This is illustrated in Fig[3jD with g — 1 and K = 1.2. In this regime, the initial population difference decays to 
a stationary value without undergoing any significant Rabi oscillations and the acousto-optic effect dominates the 
Bragg scattering effect. In this regime, we observed that since the acousto-optic effect is the dominating process, half 
the incident light energy is transferred to the diffracted beam in a short time of about t — 30ms. There are two 
additional parameters, the time dependent interaction U (r) and strength of the perturbation a which can influence 
the dynamics of the system. On increasing the interaction parameter by increasing the number of particles, the Bragg 
scattering mechanism becomes less efficient and the amplitude of the Rabi oscillations are reduced (as illustrated in 
FigEt), probably due to a dephasing caused by an increased detuning from the Bragg resonance condition. As seen 
from Eqn l21i the interaction term has a dynamical dispersive effect on the Bragg resonance. An increased number of 
particles also means that the wavenumber k s of the Faraday pattern decreases, i.e loss of the Faraday pattern. As a 
result the acousto-optic effect becomes less efficient and is not able to tranfer energy and momentum from the phonons 
to the light field. This also contributes to the observed decreased amplitude of the Rabi oscillations in Fig(3J; and 
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hence a decreased transfer of energy to the diffracted beam per Rabi cycle. The effect of increasing the strength of 
the perturbation is shown in Fig(3Ji. The beats in the Rabi oscillations are found to be enhanced, due to an increased 
contribution from the acousto-optic rate. 

Finally we would like to mention that coherence time which is basically determined by the amount of time the 
diffracted photons stay inside the condensate should be high. This can be achieved by keeping the incident angle 
9i large so that the incident photon after diffraction travels a longer path inside the condensate before it leaves the 
condensate. 



IV. CONCLUSION 

In conclusion, we have studied light scattering in the Bragg regime from an elongated cigar shaped BEC in the 
presence of a sound-like wave (Faraday waves) along the long axis. As in the case of the classical acousto-optic effect we 
show that the Faraday waves are able to transfer energy to the light field. We find that under appropriate conditions, 
there are two phase matching phenomena, the quantum acousto-optic effect (Bragg scattering of light by the dynamic 
matter-wave grating due to the Faraday waves) and the Bragg scattering of the BEC atoms from the optical grating 
formed in the BEC due to the interference between the incident and the diffracted laser beams. In section II, we have 
shown the system is relatively more unstable in the presence of the light field. In section III, we showed using the 
Heisenberg equation of motion for the atomic and light field amplitudes, that by varying the intensity of the incident 
light, one can switch between the acousto-optic regime and the atomic Bragg scattering regime. In particular, when 
the atom-field interaction is more that the light field decay rate, the atomic Bragg scattering dominates as seen by the 
Rabi oscillations between the two BEC momentum modes and in the opposite limit the acousto-optic effect dominates. 
On increasing the particle number, both the atomic Bragg scattering and the acousto-optic effect become less efficient 
as observed from the reduced Rabi oscillation and reduced amount of power transferred to the diffracted beam per 
Rabi cycle. 
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